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ABSTRACT 

We present a model for the UV to submillimeter emission from stars 
embedded in the general interstellar dust in galaxies (the 'infrared cir- 
rus' component). Such emission is characterized by relatively low optical 
depths of dust and by cool (< 30 K) dust temperatures. The model incor- 
porates the stellar population synthesis model of Bruzual & Chariot and 
the dust model of Siebenmorgen & Kriigel which includes the effects of 
small grains/PAHs. We apply the model to fit the optical to submillimeter 
spectral energy distributions (SEDs) of nearby galaxies which are domi- 
nated by cirrus emission and find that our simple model is quite adequate 
to explain the observed SEDs. 

We also, more controversially, apply this cirrus model to the SEDs of 
high redshift sources detected in blank field submillimeter surveys with 
SCUBA. Surprisingly, an excellent fit is found for many of these sources, 
with typical values for the optical depth Ay and the surface brightness of 
the stellar radiation field ip being only a factor 2-3 higher than for nearby 
galaxies. This increase is not unreasonable given the expected evolution 
of dust optical depth in currently favoured star-formation history models. 

We conclude that the tendency to interpret the high-z SCUBA galax- 
ies as very highly obscured starbursts may be premature and that these 
galaxies may be more closely linked to optically selected high redshift 
galaxies than previously assumed. 

Key words: infrared: galaxies - galaxies: evolution - star:formation - 
galaxies: starburst - cosmology: observations 



1 INTRODUCTION 

The need for several different components to un- 
derstand the infrared and submillimetre spectral en- 
ergy distributions (SEDs) of galaxies has been rec- 
ognized for some time. Rowan- Robinson and Craw- 
ford (1989) modelled IRAS galaxy colour-colour dia- 
grams in terms of three components, starburst, cirrus 
and AGN components. The cirrus component corre- 
sponded to the absorption of the general starlight of 
a galaxy by interstellar dust and reemission in the in- 
frared, as identified in our Galaxy by Low et al (1984). 
Detailed models for the cirrus component were given 
by Rowan- Robinson (1992), but these did not include 
a treatment of the line emission from the PAH com- 
ponent. A model of the infrared emission of the solar 



neighbourhood that included the effect of small grains 
and PAHs was presented by Siebnmorgen & Kriigel 
(1992). Efstathiou, Rowan- Robinson & Siebenmorgen 
(2000; hereafter ERS2000) have given the most de- 
tailed model to date of the evolution of the infrared 
emission from a starburst. In this paper we describe a 
tool which allows the cirrus emission to be modelled 
from 0.1-1000 ^m and apply it both to local galaxies 
which have been mapped with SCUBA and to high 
redshift galaxies found in deep SCUBA surveys at 
850 fj,m. Hitherto there has been a tendency to inter- 
pret these SCUBA galaxies as high redshift versions 
of Arp220, but Rowan-Robinson (2001) has pointed 
out that submillimetre counts and background radia- 
tion are much more naturally understood in terms of 
cirrus-like emission. 
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Silva et al (1998) and Granato et al (2000) have 
given models for infrared SEDs which include both 
starburst and cirrus components. These models take 
into account the distribution of dust and stars in the 
galaxy and incorporate a population synthesis model. 
The effect of changes in metallicity during the lifetime 
of a galaxy is also taken into account. The cirrus mod- 
els we present here are simpler than those of Silva et 
al (1998) and Granato et al (2000) but appear to be 
just as good a representation of the observed SEDs of 
nearby normal galEixies. Starburst models that take 
into account the distribution of dust and stars were 
also presented by Kriigcl & Sicbenmorgcn (1994) and 
Siebenmorgen, Kriigcl & Laureijs (2001). 

The layout of this paper is as follows: section 
2 describes our cirrus model in detail, section 3 ap- 
plies this model to nearby galaxies which have been 
mapped with SCUBA, section 4 discusses the appli- 
cation to high redshift SCUBA galaxies, and section 
5 gives our discussion and conclusions. 



2 THE CIRRUS MODEL 

The ingredients of a cirrus model are (1) the input 
stellar radiation field, which has to be characterized 
both with regard to its spectrum and its intensity, (2) 
an interstellar dust model, with assumptions about its 
distribution and opacity, (3) a radiative transfer treat- 
ment of the interaction between the two to generate 
the output SED. Our goal is to provide a versatile 
tool which can be used for example in simulations of 
the galaxy formation history. 

For specifying our input stellar radiation field we 
have used the Galaxy Isochrone Synthesis Spectral 
Evolution Library (GISSEL, Bruzual & Chariot 1993) 
that gives the ultraviolet to near infrared spectrum of 
radiation from a mass of stars as a function of time. 

Stars form in dusty molecular clouds and spend a 
considerable amount of time inside them. During this 
phase, the stellar radiation is absorbed by molecular 
cloud dust and reprocessed to the infrared. We use 
the code of ERS2000 to compute the emission of the 
stellar population during this phase. The molecular 
cloud is assumed to disperse 7.2 x lO'^ years after the 
instantaneous formation of stars inside it. However, 
well before 7.2 x lO'^ years non-spherical evolution of 
the molecular cloud may allow a fraction / of the 
starlight to escape without any dust absorption. We 
assume that this occurs at time tm after star forma- 
tion. By further assuming that the radiation field in a 
galax;y is due to a large number of randomly oriented 
molecular clouds, their average emission (for stars in 
the age range tm to 7.2 x 10^ years) is approximately 

(1 - f)s: + fs: 

where is the emission from a spherical CMC (as 
computed by ERS2000) and St is the emission from 



the stellar population in the absence of any molecu- 
lar cloud dust (Bruzual & Chariot). The emission of 
stars younger than tm or older than 7.2 x lO'^ years is 
assumed to be Sf, or respectively. In fact for most 
of this paper we have simplified the model by assum- 
ing / = 1. The effect of allowing / < 1 is discussed in 
section 5. 

Other parameters which have to be specified are 
the star formation history function 4>,{t), the epoch 
at which this history begins, ti, the epoch of observa- 
tion, t, and the usual parameters for the stellar IMF. 
Since the quantity which determines the dust temper- 
ature is the intensity of the stellar radiation field, we 
have chosen to characterise this in terms of the ratio 
of the bolometric intensity of the radiation field to 
that in the stellar radiation field in the solar neigh- 
bourhood (Mathis et al 1983), ijj. The star-formation 
rate, 4>»{t), is assumed to have an exponential form, 
with exponential time scale, r, and a Salpeter IMF 
from 0.1-125 Mq is assumed. The age of the galaxy, 
t* = t — ti, can also be specified. 

We use the dust grain model of Siebenmorgen 
and Kriigel (1992), which incorporates a detailed 
treatment of small grains and PAHs, and then charac- 
terize the opacity of the interstellar dust by Ay- Ay 
determines how much of the uv to near-ir light is ab- 
sorbed by dust and reemitted in the infrared and sub- 
millimetre bands, and therefore controls the ratio of 
luminosity in the far infrared to the luminosity in the 
uv to near-ir. If Ay << 1, we are in the optically thin 
regime and the dust temperature will then be deter- 
mined by ^ and the shape of the input spectrum (de- 
termined by the star formation history). For Ay > 1 
but not >> 1, the dust will modify the emergent op- 
tical and uv spectrum significantly, but will still not 
result in any significant change to the dust temper- 
ature along the line of sight to the central plane of 
the galaxy. We make a further simplifying assump- 
tion, that for the purposes of modelling the galaxy's 
emergent SED, we can characterize the galaxy by a 
single average value of Strictly speaking we should 
model the density distribution of both dust and stars 
and carry out a full radiative transfer calculation to 
compute the emergent spectrum. The emergent spec- 
trum would essentially be a superposition of spectra 
corresponding to different values of t/) at different lo- 
cations through the galaxy. In practice we find that a 
single value of tp characterizes the emergent spectra 
very adequately. 

To summarize, the parameters of the model that 
have to be set are tp,T,tt,,Ay,f and tm- The param- 
eters T,tf,f and tm determine the shape of the op- 
tical SED, with Ay determining the reddening to be 
applied to this and also, in a self-consistent manner, 
the bolometric luminosity of the infrared emission. 
ip determines the temperature of the different grain 
species in the model and hence the shape of the far 
infrared SED. There is also a slight dependence of 
grain temperature on the shape of the optical SED. 



© 0000 RAS, MNRAS 000, 000-000 



3 

UGC903 

-41. 



NGC958 

I 4 -41. I I I 




-10 



I 

NGC5962 



1 

log (V**") 
NGC6181 



§• -7 

y -8 



E 

§• -7 
? -8 



-10 



109 (X/M<ti) 

UGC5376 



log (X/>an) 
NGC2990 




-lOE 



-lOE 



log (Vaw) 



log (V»wn) 



Figure 1. Spectral energy distributions of nearby cirrus galaxies. Data from IRAS, Dunne & Eales (2001) and the 
NASA/IPAC Extragalactic Database (NED). Model parameters given in Table 1. 



For tt = 0.25 Gyrs and r = 6 Gyrs the temperature 
of the large classical grains is in the range 17-24K 
for V = 1 and 28-41K for ip = 21. Note that the 
far-infrared emission is dominated by the emission of 
the largest and coolest grains. For = 12.5 Gyrs and 
r = 6 Gyrs the corresponding temperature ranges are 
16-21K for V = 1 and 26-36K for ip = 21. 



Our code and a selection of models is available 
at http://astro.ic.ac.uk/~ane/cirrus-models.html 
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Name 


Av 


r/Gyrs 






UGC903 


0.9 


6 


4 


1.05 


NGC958 


0.4 


9 


2 


0.56 


NGC5962 


0.4 


5 


3 


0.56 


NGC6181 


0.4 


11 


3 


0.97 


UGC5376 


0.9 


5 


5 


0.35 


NGC2990 


0.4 


11 


3 


1.33 


NGC1667 


0.5 


11 


8 


1.42 



Table 1. Fitted parameters for the nearby cirrus galaxies: dust optical depth, Ay; exponential time-scale for star formation, 
t; intensity of interstellar radiation field relative to solar neighbourhood, V; goodness of fit, x^. All models assume a galaxy 
age tt of 12.5Gyrs, / = 1 and tm = 3Myrs. 



3 APPLICATION TO NEARBY 
GALAXIES 

Dunne & Eales (2001) have given the results of a 
SCUBA mapping program of 19 galaxies. The far in- 
frared spectrum of four of these (N1614, I Zw 107, 
IR 1525-1-36, Arp 220) arc obviously dominated by a 
starburst component. For a further 7 galaxies there 
is evidence from the 60/850 //m colours that a star- 
burst component is present ('opio^^^jj > 1.8). We 
also excluded NGC7541 because of the large flux cor- 
rection applied to the 450/itm data by Dunne & Eales 
(2001). We have therefore modelled the SEDs of the 
remaining 7 galaxies, including optical and near in- 
frared data from the literature, in terms of a cirrus 
component (Figure 1). The parameters for the mod- 
els arc given in Table 1. For NGC1667 we only give the 
fitted parameters in Table 1. The fits are extremely 
good and do not show evidence that a more elabo- 
rate cirrus model is required. The values of i/> range 
from 2 to 8, and the values of Ay range from 0.4 to 
0.9. Such values are consistent with the local distri- 
bution function for the ratio logw{L fir / Lopt) , which 
peaks at -0.7 and can be fitted with a Gaussian with 
G = 0.24 about this (Rowan- Robinson et al 1987). As 
pointed out by Rowan-Robinson (1992), the parame- 
ter V' is proportional to the mean surface brightness 
of the galaxy. 



4 APPLICATION TO HIGH REDSHIFT 
SCUBA GALAXIES 

The detection of a number of galaxies in blank field 
surveys with SCUBA et 850 ^.m (Hughes et al 1998, 
Eales et al 1999, Bargcr ct al 1999, etc) and the re- 
alization that most of these are at redshift > 1, with 
some probably at z > 3, has had a big impact on ideas 
about galaxy formation and evolution. Most attempts 
to fit the SEDs of these galaxies (Hughes et al 1998, 
Downos et al 1999, Lutz ot al 2001, etc) and most 
attempts to fit the 850 yum counts and background 
(Guiderdoni et al 1998, Dole et al 2001, Xu et al 
2001, Franceschini et al 2001, Pearson 2001, Granato 
et al 2001) have assumed that we are looking at a 
very strongly evolving, deeply dust-embedded star- 



burst population. However Rowan-Robinson (2001) 
argues that the natural interpretation of the submil- 
limetre counts and background is in terms of a cirrus- 
type component. How can the latter picture be recon- 
ciled with the failure to identify many of the SCUBA 
galaxies with optical counterparts ? Obviously the 
dust optical depth must be higher than in local disk 
galaxies. However that is not an unreasonable expec- 
tation a.t z — 1-3. Calzetti and Heckman (1999) and 
Pei et al (1999) have predicted how the typical dust 
opacity in galaxies should evolve with epoch for star- 
formation histories typical of those derived from uv 
and infrared surveys. An increase of the average Ay 
in the interstellar medium in galaxies by a factor of 
2 or 3 over present-day values at z =1-2 is to be ex- 
pected, because the increased gas-density in galaxies 
far outweighs the slight decline in metallicity over the 
same look-back time. 

Here we have selected two samples of SCUBA 
high-z galaxies from SCUBA blank-field surveys, con- 
fining attention to those with the most reliable optical 
identification: (a) submm galaxies which have been 
confirmed by submillimctrc intcrforomctry (e.g. Lutz 
et al 2001), (b) submm galaxies from the 8 mjy sur- 
vey (Scott et al 2001, Fox et al 2001) which have been 
confirmed by association with radio sources (Ivison et 
al 2002). For both these samples the optical and near 
infrared associations should be reasonably secure. For 
sample (a) we have the added bonus that the reality 
of the submillimetre source is confirmed. 

To make use of the radio data we have ex- 
tended our model to the radio regime by utilizing the 
well known far-infrared-radio correlation (Helou et al 
1985). For A > 60/im we add to the rest frame model 
SED S:, where 

5^^ = 2.7x10 ^{2.58S6o^m + Sioonm)x{ ^ AGHz ^ 

where 560m"i and S'ioomwi are the 60 and 100/im fluxes 
predicted by the model and i/ is the frequency in Hz. 
The factor 2.7 x 10"^ is the mciximum value allowed 
by the far- infrared-radio correlation (Blain 1999). 

The model allows us to self-consistently estimate 
the star formation rate (SFR) of the objects studied. 
Assuming f« = 0.25Gyrs the SFR is given by 
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Figure 2. Spectral energy distributions of high redshift cirrus galeixies confirmed by millimetre interferometry. Data from 
Lutz et al (2001), Dunlop et al (2002), Ivison et al (1998, 2000, 2001), Prayer et al (2000) and Gear et al (2000). Model 
parameters given in Table 2. For HDF850.1 only the z = 4.5 fit is shown. 



7.88 X 109 ^ 

where Lboi is the bolometric luminosity of the ob- 
ject (in units of solar luminosity) which we obtain by 



integrating over the whole wavelength range covered 
by the models. 

To test the predictions of our cirrus model we 
generated a grid of models in which we vary ip (in 
the range 1-21 in steps of 1) and Ay (in the range 
to 3 in steps of 0.1). For most of the sources in these 
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Figure 3. Spectral energy distributions of high redshift cirrus galaxies from the 8 mjy SCUBA survey in ELAIS-N2. Data 
from Fox et al(2001), Ivison et al (2002). Model parameters given in Table 2. 



samples the redshift is unknown so we treat it as a free 
parameter varying between and 10 in steps of 0.1. 
For objects with data in only four bands (all of the 

objects in Lockman East and SMM00266+170) we fix 
tp to 6. For all the models we assume t« = 0.25Gyrs, 
T = 6 Gyrs, / = 1 and tm = 0. 

Minimum fits of our cirrus model to the SEDs 



of these two samples of high z SCUBA galaxies are 
shown in Figures 2-4, with the parameters of the mod- 
els given in Table 2. In N2 we model all the objects 
with robust radio associations except N2850.5 which 
is a blank field in all optical and near-IR bands, and 
N2850.13 which is detected only in K. In LE we ex- 
clude LE850.12 and LE850.18 which are detected only 
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Figure 4. Spectral energy distributions of high redshift cirrus galaxies from the 8 mjy SCUBA survey in Lockman. Data 
from Fox et al(2001), Ivison et al (2002). Model parameters given in Table 2. Fits with a cirrus model that includes a dust 
enshrouded phase, discussed in section 5, are shown with the dotted lines. 



in the I band. LE850.1 is included in the sample with 
millimetre interferometry. Fits that exceed the upper 

limits are rejected. The quoted errors in z indicate 
the spread in z for fits with less than Xmin + 1- 

The fits we obtain are generally good. One no- 
table exception is SMM02399-0136 where the high 



reduced is almost entirely due to the radio point. 
Ivison et al. (1998) find evidence for an AGN in this 

object. Also, the high fitted value of i/) in this ob- 
ject (the max;imum value in our grid of models) is 
perhaps indicative of the presence of a circumnuclear 
starburst accompanying the AGN in this object. The 
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7.9; 
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0.8 
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1.541 


(6) 


0.7 
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0.086 


(6) 


0.2 
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(6) 
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s.i; 
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78.395 


(6) 


0.8 
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lensing logio 
magn. 



+0.5 
-0.5 
+0.1 
-0.1 
+0.1 
-0.1 
+0.3 
-0.2 
+0.2 
-0.1 
+0.1 
-0.1 
+0.1 
-0.2 



+0.1 
-0.1 
+0.1 
-0.1 
+0.1 
-0.1 
+0.1 
-0.1 
+0.1 
-0.1 
+0.1 
-0.1 



+0.4 
-4.3 
+1.1 
-0.6 
+0.2 
-0.1 
+0.4 
-0.2 
+0.3 
-0.1 
+0.1 
-0.1 
+0.1 
-0.1 



+0.2 
-0.4 
+0.1 
-0.1 
+0.3 
-0.2 
+0.1 
-0.5 
+0.3 
-0.3 
+0.1 
-0.1 
+0.1 
-0.1 



3 

2.5 

3 
3 

2.4 



12.76 
13.31 
13.05 
12.23 
12..53 
12.81 
12.91 



12.45 
12.75 
13.22 
13.01 
12.74 
12.63 



13.46 
12.72 
13.02 
12.81 
13.09 
12.69 
13.27 



13.52 
13.11 
13.08 
12.95 
13.16 
12.72 
13.30 



SFR 



731 
2577 
1441 
219 
436 
819 
1033 



360 
705 
2121 
1295 
700 
542 



3625 
673 

1343 
815 

1566 
620 

2366 



4222 
1627 
1522 
1141 
1822 
662 
2545 



Table 2. Fitted parameters for the high redshift cirrus galaxies with fixed values given in brackets. The age of the starburst 
t* is fixed at 0.25 Gyr and the star formation history is assumed to be exponentially decaying with time constant r = 6 
Gyrs. Unless otherwise stated we assume / = 1 and tm = 0. Luminosities are calculated by integrating over the whole 
wavelength range covered by the models. We assume a flat Universe with Ho = 65 Km/s/Mpc and A = 0.7. Where 
appropriate the luminosities (and the resulting SFR) have been corrected by the lensing magnification factor given in the 
table. 



fits to LE850.3 and LE850.7 (not plotted) are also 
poor. The steep radio spectrum of tliese objects may 
also suggest tlie presence of an AGN (Ivison et al 
2002). 

Tiie values of tp arc higher, typically, by a factor 
of 2-3 than the values for local galaxies, and the values 
of Av are also higher by a factor of 2-3. Although the 
values of Ay are higher than in local galaxies, they are 
much lower (by one or two orders of magnitude) than 
would be expected in a typical molecular cloud under- 



going a starburst. So these models are radically differ- 
ent from the typical dusty starburst model for these 
high-redshift submillimetre galaxies. Of course, be- 
cause the bolomctric luminosities correspond to high 
rates of star-formation, much of the optical and uv 
light which is illuminating the interstellar dust in our 
cirrus models is light that has escaped from a star- 
burst region, for example by non-spherically symmet- 
ric evolution of the associated HII regions. It is also 
interesting to note that the Ay we infer are compa- 
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rable to those of molecular clouds in late stages of 
their evolution (about a few times 10^ years) as pre- 
dicted by the starburst evolution model of ERS2000. 
As pointed out by the latter authors the predicted 
IRAS colours of these clouds are similar to those of 
cirrus galaxies. 

With a few exceptions we are able to put fairly 
strong constraints to the redshifts of the sources we 
model. For SMMS14011+0252 Ivison et al. (2000) re- 
port a spectroscopic redshift of 2.56 which is consis- 
tent with our model. Small et al. (2003) give a spec- 
troscopic redshift of 2.38 for N2850.4 which is close 
to the deduced photometric redshift. For HDF850.1 
there are two distinct minima, one at z = 4.5 and 
one at around z = 7.1. In Table 2 we give the fitted 
parameters for both models. For LE850.2 there is a 
broad minimum between 2 = 4 and 9. 

Of the 23 objects in our samples, our models suc- 
cessfully fit the SEDs of 16 of them (70 %). For 4 
objects (17%) we have insufficient data (fewer than 
4 bands) to be able to test our model. For these ob- 
jects we may need a higher optical depth model of the 
M82 or Arp 220 starburst type or values of Av (> 3) 
from our cirrus model which are not really consistent 
with the ^ Ay not >> 1' requirement of our model. 
For 3 objects our fitted i^) is close to the maximum in 
the range considered (21). We assumed the age of the 
starburst = 0.25 Gyr in each case. This will under- 
estimate the contribution of older, low-mass stars in 
the rest-frame near infrared if present. In the case of 
N2850.8 we do see evidence for such a near-IR excess. 



5 DISCUSSION AND CONCLUSIONS 

Although the current state of observations of high- 
z SCUBA galaxies is insufficient to tell whether our 
models are better or worse than starburst models, the 
fits shown here demonstrate that it is premature to 
make an analogy between SCUBA galaxies and dusty 
ultraluminous starburst galaxies like Arp 220. Quite 
a modest dust optical depth [^Ay ~ 1 — 3) combined 
with the redshifting of the optical and uv radiation 
is sufficient to explain why these galaxies are faint in 
the optical. 

The inferred bolomctric luminosities imply that 
enhanced star formation has taken place in these 
galaxies, so even if we are correct that the observed 
submillimotro omission arises mainly from cirrus emis- 
sion, ie reemission from dust in the general interstel- 
lar medium of the galaxy, it is likely that some op- 
tical and ultraviolet light is intercepted by dust in 
the dense molecular gas from which the stars formed. 
Such emission would peak at ~ 50 (1+z) /im and 
could eEisily contribute bolometric luminosities com- 
parable to the cirrus luminosities of table 2 without 
being detectable at present. Estimates of the fraction 
of optical-uv light which escapes from star-forming re- 
gions in the local universe range from 25-50 %. In Fig 



4 we have also shown fits for a more elaborate model 
in which newly formed stars arc completely shrouded 
in their molecular cloud for the first lO'^ yrs, 50% of 
the starlight escapes to illuminate the interstellar dust 
from 10'^ - 7.2 x 10** yrs (i.e. / = 0.5 and = 10^ 
years) and 100% escapes thereafter (cf the starburst 
evolution models of ERS2000). The starlight that es- 
capes is reprocessed to the infrared in the same way 
as the pure cirrus model. We have generated a grid 
of models spanning the same parameter space in Ay, 
ip and z as for the pure cirrus model. These mod- 
els, when observed at 0.25 Gyr from the start of the 
starburst, arc very similar to our pure cirrus models 
except for additional radiation peaking at 50 /im. 

A critical test of our models would be to mea- 
sure the angular extent of the millimetre or submil- 
limetre sources in high-z galaxies. In our model these 
should bo of order an arcsoc, within range of cur- 
rent millimetre-wave interferometers, whereas the ob- 
scured nuclear starburst model would imply sizes an 
order of a magnitude smaller. Of course, this test 
would not be able to discriminate between our model 
and a model whore star formation takes place in an 
extended highly obscured starburst. 

There is a further factor at work which exagger- 
ates the difference between the SCUBA galaxies and 
the high-redshift galaxies found in Lyman break sur- 
veys (Steidel et al 1999). Assuming that there is a 
spread in the typical optical depth of the interstellar 
medium, uv surveys will inevitably be biased towards 
galaxies with lower than average optical depth, while 
submm galaxies will be biased towards higher optical 
depth. According to this picture, if the Lyman break 
surveys could go a factor 10 deeper they would start 
to see the SCUBA galaxies, and if the SCUBA sur- 
veys could go a factor 10 deeper they would start to 
see the Lyman break galaxies. This is already hinted 
at by the analysis of Peacock et al (2000) . 
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